Fast approximation for the energy and density of electrons in two dimensions 
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We derive a self-consistent variant of the energy functional presented in S. Pittalis and E. Rasanen, 
Phys. Rev. B 80, 165112 (2009). The resulting scheme is an exphcit orbital-free representation 
of the electron density and energy in terms of the external potential, the number of electrons, and 
the chemical potential determined upon normalization. The method yields reasonable results for 
two-dimensional quantum dots, rings, and slabs with a practically negligible computational cost. 
Both the benefits and drawbacks of the approach are analyzed. 
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Orbital-free density-functional theory (OF-DFT) is a 
computationally appealing method to deal with large sys- 
tems beyond the reach of conventional DFT. At present, 
OF-DFT methods allow to deal with systems up to a 
million atomsi These methods await to be fully ex- 
plored in the context of low-dimensional systems and na- 
noelectronic devices. In two-dimensional (2D) physics 
one of the main challenges of DFT is to deal with re- 
gions of the 2D electron gas^ comprising hundreds or 
thousands of interacting electrons, e.g., in the quan- 
tum Hall regime. During the recent years successful ef- 
forts have been made to describe quantum Hall devices 
by incorporating the electron-electron (e-e) interactions 
up to the TF level, but improvements regarding, e.g., 
the exchange-correlation effects neglected in the TF ap- 
proach are needed. 

As the name suggests, OF-DFT is free from the use of 
the Kohn-Sham orbitals needed in the calculation of the 
Kohn-Sham kinetic energy, and thus the only explicitly 
needed variable is the electron density p{v). The earliest 
OF-DFT method can be dated back to the TF theory 
employing the exact result of the homogeneous electron 
gas for the kinetic energy, and the Hartree approximation 
for the e-e interaction. In fact, most orbital- free schemes 
can be regarded as modifications or improvements to the 
TF method.i. The crux of the problem in OF-DFT is to 
find an approximation for the (non-interacting) kinetic- 
energy functional that may be generally applicable. For 
this ambitious goal, a promising approach is an orbital- 
free formulation exploiting the potential rather than the 
the density as the basic variable. - As described below, 
in this work our aim is not to improve over the TF ap- 
proximation for the kinetic energy.^ Instead, our focus 
is on replacing the Hartree energy with a much simpler 
expression accounting also for quantum part of the e-e 
interaction to some extent. 

Recently, a 2D orbital-free expression for the energy 
was shown to lead to a major improvement over the TF 
results when applied to quantum dots and slabs up to 
200 interacting electrons.* However, the energies were 
not obtained self-consistently, but by using the electron 
densities from the local-density approximation (LDA). 



In this Brief Report, we show that the approach works 
well also when it is applied fully self-consistently. In do- 
ing so, we express the density in terms of the given ex- 
ternal potential, number of electrons, and the chemical 
potential - the latter determined through the particle 
normalization condition. The functional is tested for a 
variety of 2D systems including quantum dots, rings, and 
slabs. A promising overall performance is obtained with 
an almost negligible computational cost. This gives the 
method an appealing position in the field of OF-DFT for 
low-dimensional systems and electronic devices. 

In Ref. H we showed that using the Gaussian approxi- 
mation for the cylindrical average of the pair density, and 
a uniform coordinate scaling, the total e-e interaction en- 
ergy can be expressed as 

WW)] = / drp3/2(r), (1) 

where N is the number of electrons and p{y) is the 2D 
total electron density. In turn, an orbital-free estimation 
for the total energy can be expressed as 

Etot[p{r)]=TTF[p{v)]+W[p{v)]+ j drp(r)vext(r) (2) 

where 

TtfW)] = \j dvp'iv) (3) 

is the TF kinetic energy in 2D, and the last term is the 
energy contribution due to the external scalar (confining) 
potential Wcxt(r). 

Now, we emphasize a few important points of the en- 
ergy functional given in Eq. First, for TV = 1 the 
functional evidently reduces to the noninteracting TF ap- 
proximation. Second, for N ^ 1 the interaction contri- 
bution is a single term that has the simple form of the 
exchange energy in the LDA but with a different prefac- 
tor y/N - 1. In this way, the computational burden of 
the Hartree term in the TF approximation has been elim- 
inated. On the other hand, W[p] in Eo;^ ((!]) is evidently 
not size-consistent as discussed in Ref. [3. 
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In order to find the ground-state density we have to 
minimize Eq. ^ for a fixed number of particles. We 
may first express the total energy in a single integral as 



^tot[p(r)] = I drF[p{v)] ^Jdr 



2P 



TT N -1 



(4) 



We have to find a stationary value for the functional 
i^[/o(r)] with respect to variations in p{r). To take the 
electron number conservation into account, we introduce 
another functional G'[p(r)] — p{r) so that 



J drG[p{r)]=N. 



(5) 



This constraint introduces a Lagrange multiplier p in the 
variational equation, which can be written as 



dF _^ 
dp ^ dp 

Substituting F and G to this equation yields 



7rp(r) H ,, 

^ 4 V 2 



(6) 



(7) 



As this expression is quadratic in p^^^, we find an explicit 
expression for the density, 



p(r) 




N -I 
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OZ TT 



(8) 



This is our key result showing that the density can be 
solved instantaneously for any external potential Ucxt and 
any N. The only variable to be determined numeri- 
cally is p that follows from the normalization condition 
in Eq. ([5]). The symbol [...]^ in Eq. ([8]) represents an ad- 
ditional constraint that no sign changes under the square 
in Eq. ([8]) (leading to unphysical "nodal lines" in the den- 
sity), nor negative values under the square- root (leading 
to complex densities), are allowed. Once p{r) is deter- 
mined through Eq. the total energy is obtained from 
Eq. ©. 

Let us emphasize the difference between the present 
and and the TF approximation. In the latter, the varia- 
tional procedure applied to the total energy leads to an 
integral equation for the density. The TF scheme then 
transforms into a differential equation (which in 3D leads 
to the Poisson equation). Instead, our functional is free 
from this complexity due to the simple expression for 
the interaction energy [Eq. ([Ij] in comparison with the 
Hartree integral utilized by the TF method. Although 
the Hartree term is simple to calculate in most appli- 
cations, it may become a bottleneck in very large 2D 
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FIG. 1: (color online) (a) Cross section of the electron density 
in a square quantum well with side length L = 1 at j/ = 0.5. 
The well contains A'^ noninteracting electrons, (b) Relative 
error in the total energy given by the present functional for 
the same system as a function of A'^. 



systems. On the other hand, the numerical cost of the 
present scheme is practically negligible for any N . 

Before computational applications we consider some 
analytic limiting cases. First, as noted above, the in- 
teraction energy correctly vanishes at A'^ = 1. On 
the other hand, the density in Eq. (|8]) becomes simply 
PAr=i(r) = [p — Vexti^)]+/'!^, i-S., the density profile is the 
additive inverse of the external potential, and the height 
of the profile is determined by p such that the normal- 
ization condition is satisfied. As a simple example, for a 
2D harmonic oscillator the exact Gaussian-type density 
is approximated by a parabolic one. 

Second, let us consider a 2D square quantum well with 
side length L = 1 and infinite boundaries consisting of A^ 
noninteracting electrons. Thus, we omit W in Eq. ^ and 
apply the remaining expression self-consistently. The ex- 
act density can be written simply as a sum of a (squared) 
product of sine (or cosine) functions. In Fig. [ija) the 
cross section of the exact density per electron at y = 1/2 
is compared with the present functional that yields a 
constant density p = N, a.s t'oxt(r) = inside the well. 
Clearly, our approximation for the total density (dashed 
line with p/N = 1) improves as a function of A^. This 
favorable property is refiected in the estimation of the 
total energy: as shown in Fig. [IJb), the energy - here 
corresponding solely to the TF kinetic energy - becomes 
more accurate as a function of A^ in comparison with the 
exact result, i.e., the sum of the energy eigenvalues. 

Next we test the functional for a set of 2D quantum 
dots, rings, and slabs including the e-e interactions. We 
use the 2D LDA^ in DFT as our reference method; in 
the range of systems and parameters considered here the 
LDA has been shown to provide ~ for the present pur- 
pose - sufficiently accurate total energies and densities 
(see, e.g., Refs. ^^ and ITTI for quantum dots and rings, 
respectively). The LDA calculations are performed using 
the octopus codej^ 

Figure [2ja) shows the radial densities obtained from 
Eq. ^ for 2D harmonic quantum dots defined by 
Vcxtir) = Ld'^r'^/2 with ui = 0.5 (solid lines). The dashed 
lines show the corresponding LDA results. As expected. 
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FIG. 2: (color online) (a) Electron densities in harmonic quan- 
tum dots containing A*' = 6 and 56 electrons, respectively. 
The dashed lines denote the conventional density-functional 
calculations with the local-density approximation. The solid 
lines correspond to the results of the present orbital-free func- 
tional, (b) The same as in (a) but for two quantum rings 
containing = 12 and 38 electrons, respectively. 



the present functional cannot describe the shell structure 
due to the lack of orbitals. However, the overall shape of 
the density profile is found in a satisfactory fashion, and, 
most importantly, the correct radial extent of the density 
profile is obtained in both cases {N = 6 and N = 56). 
Note that the correct density has a Gaussian decay^^'^ 
whereas in our approximation, by definition, the decay 
follows the external potential up to the power of one. As 
an obvious drawback, agreement with the LDA does not 
seem to improve as a function of N. 

In Fig. mja) we show the corresponding result for 
a quantum ring modeled by Ucxt(''') = -|- 
Vbexp(-rVd2) with w = 0.5, Vq = 20, and d = 1. The 
model potential is the same as the one used in Refs. |n|, 
[13, and [H, the last reference showing direct comparison 
with an experiment. We find a reasonable qualitative 
agreement between the present functional and the LDA. 
The qualitative agreement is similar for both = 12 and 
A^ = 38. 

Apart from densities, it is important to assess the per- 
formance of the present functional for total energies. In 
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FIG. 3: (color online) (a) Relative error in the total energy 
(with respect to the local-density approximation) for rectan- 
gular quantum slabs with = 12 . . . 200 as a function of the 
density parameter Va. The results of the present functional 
and that of Ref. @ - employing a non-self-consistent version of 
the present functional - have been indicated, (b) The same 
as (a) but for the electron-electron interaction energy. 



Ref. 8 this was done non-self-consistently by applying 
the LDA densities as an input for the total-energy func- 
tional in Eq. ©. Here we use the present self-consistent 
approach to the same extensive set^ of rectangular quan- 
tum slabs li^iii The external potential is defined by steep 
boundaries, so that the side-length ratio is fixed to two, 
and the shorter side is varied between 5 and 60 a.u. The 
number of electrons is varied in the range A^ = 12 . . . 200. 

In Fig. ^a) we show the relative total-energy dif- 
ferences from the reference LDA results for a set of 
rectangular quantum slabs as a function of the average 
electron density in terms of the Wigner-Seitz radius?^ 
Ts = \/ A/ (ttN), where A is the area of the slab. Results 
are shown for the present functional and for the non-self- 
consistent calculation of Ref. H. The latter scheme seems 
to outperform the former at r<, > 4 when the average 
density becomes low. This improvement is accidental in 
the sense that the self-consistent density of our functional 
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[Eq. ([8])] remains constant in the slab regardless of A'' and 
Tg, as Vcxt is zero inside the boundaries. As a function 
of N, however, it can be expected that the present func- 
tional becomes more accurate, since the actual density 
becomes flatter in the large-TV limit - similarly to the 
noninteracting case analyzed within Fig. [1] above. 

Figure 131b) shows the relative errors in the interaction 
energy W in Eq. ([Ij that can be compared to the LDA re- 
sult, Wlda = Eh + Exc, where xc refers to the exchange 
and correlation. The results are similar to those for the 
total energy, though now the present self-consistent func- 
tional seems to perform slightly better than the non-self- 
consistent one. 

Quantum slabs as defined above are also useful model 
systems in the analytic and self-consistent application 
of our approximation in limiting situations. In contrast 
with the analysis above (Fig. [1]) we now include the e-e 
interactions. First, when the area A of the slab is fixed 
and TV goes to infinity we get ^ [Tr/A + 3/(4^2)] N 
and -Btot ~ (a*/2 + N/8)N. Then, when iV >> 1 but 
finite and A goes to infinity, we get ^ w 3nN/{4:^/2A) 
and E'tot ~ 2/iA^/3. Finally, when both N and A ap- 
proach infinity while keeping the (self-consistent) den- 
sity constant, we get, again, ^ « Btt A^/(4-\/2]4) and 
Etot ~ 2/iiV/3. In all these three limiting cases, our ap- 
proximation apparently fails to satisfy the fundamental 
relation /i = dEtot/dN. We attribute this fact to the fact 



that W is not size-consistentiS. In our preliminary tests 
we have not found severe consequences of this deficiency: 
in quantum dots, for example, the relative errors in en- 
ergy differences E{N) — E{N — 1), when compared to the 
LDA calculations, are of the same order as the relative 
errors in the total energies themselves. 

Despite the obvious formal deficiencies, the present ap- 
proach has a few important advantages. First, the ap- 
proximation is extremely fast in the sense that no itera- 
tions apart from the numerically trivial determination of 
fi are required. Hence, the functional could be used as the 
first approximation for the electron density in the con- 
ventional Kohn-Sham scheme. Secondly, despite the re- 
markable simplicity the functional yields reasonable mag- 
nitudes for the total and interaction energies. Hence, the 
present analysis supports further improvements and ex- 
tensions in developing orbital-free methods for large-scale 
two-dimensional electronic devices. 
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